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INFRARED RADIATIVE TRANSFER THROUGH 


A REGULAR ARRAY OF CUBOIDAL CLOUDS 
Harshvardhan and Janies A. Weinman 

ABSTRACT 

A study has been made of infrared radiative transfer through a regular array of cuboidai clouds 
wliich considers the interaction of the sides of the clouds with each other and the ground. Tlie theory 
is developed for black clouds and is extended to scattering clouds using a variable azimuth two-stream 
(VATS) approximation (Harshvardhan et ai., 1981). It is shown that geometrical considerations of- 
; ten dominate over the microphysical aspects of radiative transfer through the clouds. For example, 
the difference in simulated IQ/sm brightness temperature between black isothermal cubic clouds and 

i 

< cubic clouds of optical depth 10, is less than 2** for zenith angles less than 50° for all cloud fractions 
when viewed parallel to the array. 

The results show that serious errors are made in flux and cooling rate computations if broken 
clouds are modeled as pianiform. Radiances computed by the usual practice of area-weigiiting cloudy 
and clear sky radiances are in error by in briglitness temperature for cubic clouds over a wide 

range of cloud fractions and zenith angies. It is also shown that the lapse rate does not markedly 
affect the exiting radiances for cuboidai clouds of unit aspect ratio and optical depth 10. 
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INFRARED RADIATIVE TRANSFER THROUGH 
A REGULAR ARRAY OF CUBOIDAL CLOUDS 


1 . Introduction 

Radiative parameterization of clouds is required for climate and general circulation models be- 
cause clouds represent the most variable quantity that affects radiative fluxes. To date, all general 
circulation models have utilized plane parallel theory to construct cloud parameterizations even 
though it has been realized that a significant portion of the total cloud field is broken and individual 
cloud elements have horizontal dimensions of the same order as their height. In these cases, the 
assumption of infinite horizontal extent in solving the radiative transfer equation is invalid and 
parameterizations based on plane parallel tlieory must be regarded as suspect. 

Tlie particular case of monochromatic radiative transfer through single cuboidal clouds has been 
studied in both solar and infrared wavelengths (McKee and Cox, 1974; Davies, 1978; Liou and Ou, 
1979 and Harshvardhan et al., 1981). However, each element in a broken cloud field cannot be 
considered independently as there is cloud-cloud interaction. The problem of infrared transfer 
through an array of black clouds lias been dealt with extensively in the Estonian literature (Ohvril 
and Epik, 1978). EUingson and Kolczynski (1980) have computed infrared heating rates in the 
atmosphere in the presence of an array of black cylindrical clouds. Interpretation of infrared radi- 
ances measured by satellite borne sensors also requires an understanding of the effects of the finite 
horizontal extent of clouds. 

The present study considers a regular array of identical cuboidal clouds, as in Aida (1977), over- 
lying a non-reflecting surface embedded in a non-participating atmosphere. Only infrared radiative 
transfer Is considered. Tlie theory of transmission through the array and emission from the cloud 
field is first developed for black isothermal clouds. The work is then extended to 1 Cl radiation 
interacting with a field of clouds which scatter and absorb radiation. Tlv? VATS solution to the 
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single cloud problem employed by Harshvardhan et al. (1981) is used. Computations to derive both 
fluxes and radiances are developed as a function of the cloud fraction. 


2. Radiative Transfer through an Array of Black Qouds 

Consider an idealized broken cloud Held in tlie form of an extended regular array of cuboids. 

A portion of this array is illustrated in Fig. 1 . Note that all elements of the array are identical and 
lie in one plane. Also, the horizontal dimensions of all clouds are assumed equal, although this con- 
straint may, in principle, be relaxed. In general, the height of each element, z*, is different from 
the horizontal dimension, s. From the figure, tlie area of the top and bottom face is Aq » s^ and of 
each side face, Aj * sz*. We define the aspect ratio of each element as 


A, 


( 1 ) 


If the spacing between each element is equal in the x- and y-directions as shown in the figure, 
then the fraction of the field covered by clouds when viewed normally, lienceforth referred to as 
the cloud fraction is 



Consider a radiatively black cuboid at uniform temperature, Tq. The monochromatic power 
(W/am"^) emitted to space by the top face is ffB(jAo where B,, is the Planck emission at T,,. The 
monochromatic power emitted by each side face is itBqAj and for a single cloud, half escapes to 
space and half goes to the ground. It is assumed that the atmosphere does not interact witii the 
radiation field. This would be reasonably true in the atmospheric window centered around 10 jam. 
The total power emitted to space by a single cuboid is 

P»7rBoAo(l +2a) (3) 

where a is the previously defined aspect ratio. We define the quantity in parenthesis as tlie power 
ratio , p. However, if the cloud is part of tlie array shown in Fig. 1, some of the energy emitted by 
the side faces will be intercepted by neighboring clouds. If the fraction of energy diffusely emitted 
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by one side face that is intercepted by all these clouds is 5, then the pC/wer emitted to space by each 
side face is jtBqAj (1 “5)/2 and the total power emitted to space by each cuboid is 

P-ffBoAo[l +2a(l -3)1 (4) 


so tliat the power ratio is 


p* 1 +2a(l -3) 


(5) 


We now introduce the concept of effective emittance for the cuboidal array on the same lines as the 
emittance in plane parallel theory. Here we use power rather than flux because the emitting areas 
have different dimensions. We define the effective emittance of a broken cloud array as the ratio of 
the power emitted by a region containing the array in the absence of ground emission to the power 
that would be emitted if the entire area was covered by a black surface at the same temperature. If 
there are n identical clouds of tops Aq as in the present model, then the effective emittance is 


nP 


. , irBoAt 

where A^ U the total area; but the cloud fraction is N » nA^/Aj and therefore from (4), (5) and (6) 
we have 


( 6 ) 


e-Np 


(7) 


Consider some limiting cases of 'the above equation. When a 0, p 1 and S * N and the problem 
reduces to the planiform case. When N 1 the clouds in the array approach each other and all the 
energy emitted by the side faces is intercepted, so 5 -*■ 1 and p -*• 1 ; here, S » N * 1 . It can be seen 
that S plays the role of an effective cloud fraction, N^, and the two can be used interchangeably for 
black isothermal clouds. 


The effective transmittance through the array of black cuboidal clouds is 

3»1-8=I-Ng (8) 

Thus, the power emitted to space from a broken cloud array at Tq over a black ground at tempera- 
ture Tj is 
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Pqui * 


( 9 ) 


where Bj is the Pbnck emission corresponding to . Further, the flux per unit area of the field is 
given by 

p 

-rrBoNe+rrBjd-Ne) (10) 

and the cuboidal black isothermal cloud problem takes the form of the planiform cloud model for 
which 

Fout’^‘BoN + ^BlO-N) (11) 

a. Effective Qoud Fraction 

If the sides of the clouds are black emitters and the clouds are isothermal, the faces will radiate 
diffusely and uniformly. The fraction of energy emitted by each side face that Is intercepted by all 
the neighboring faces can be computed using the angle factors described by Sparrow and Cess (1978). 
A top view of the regular array is shown in Fig. 2 . Due to symmetry we need only consider one side 
face, such as the hatched face on cloud 00 marked in the figure. A fraction of the energy emitted by 
this face is intercepted by the opposite face (hatched on cloud 01). All of the other hatched faces 
(as well as their mirror images) also intercept some fraction of this emitted energy. These fractions 
can be computed using angle factor tables and related equations given in the Appendix of Sparrow 
and Cess (1978), as long as the faces are not partially obscured. Even when there is some obscura* 
tion, the faces can be divided into smaller strips so that only the inradiated areas contribute to the 
total conflguration factor. This is best accomplished graphically for the different cloud spacings 
corresponding to different cloud fractions. This method of obtaining 5 and hence, N^, is prone to 
considerable error as N increases so an alternative technique was used. 

Tlie effective cloud fraction, N^, was deduced experimentally from the transmission of a diffuse 
ligiit source through an opaque array. This is the optical analog obtained by setting B(, = 0 in ( 1 0). For 
this, an array of black blocks was placed on a light table covered with a diffusing sheet. The flux 
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passing througit the array was measured with a Tektronix Model J16 digital photometer with a cosine 
corrected remote illuminance probe. The probe was translated above the blocks and the mean of a 
number of readings was used to compute the transmission, (I - N,), from (10). As a check on the 
accuracy of the experiment, the measured results for an array of infinitely long bar clouds was com- 
pared with theory which is readily available in Sparrow and Cess (1978) for this configuration. Tliis 
is shown in Fig. 3 for bar clouds of unit aspect ratio. The measurements were taken by lining up a 
large number of black blocks in a parallel array to simulate the infinite extent 'n one horizontal 
dimension. Four sets of measurements were taken and they .ire sitown in the figure with their re- 
spective uncertainties. Correspondence is good for N ^ 0.6, but for higher values, grazing li^t re- 
flected from the side walls of tlie blocks which were not perfectly black. This produced higher values 
of the measured transmission which gave an underestimate of N,. Figure 3 also shows the curve N 
corresponding to planiform clouds with a * 0. 

Measurements were then made for tlie regular array of cuboids shown in Fig. 2 with two differ- 
ent aspect ratios; the resulting values are shown in Fig. 4. For comparison, the values derived from 
summing up the angle factors is shown for small 14 when there is not much obscuration and die theor;; 
is sdll quite accurate. If the blocks are arranged such that their comers are touching in a checkerboard 
pattern, then N > 0.5 and can be obtained exactly from the expression for angle factors because 
each side face sees only three other faces. The comparison for this case is also shown on Fig. 4. 

Finally, a series of curves are empirically fitted through the available data points to represent N^, 

The empirical expression 

, [1.^2a (l ^0,15N)1N 

® 1 +2aN(! +0.15N) ^ 

is shown in Fig. 4. The curves follow the general pattern of the results presented in Ohvril and Epik 
(1978) and Ellingson and Kolczynski (1980). It may be noted that as N -*• 0, Ng (1 + 2a) N, which 
is the correct limit for a single cuboid (see Eq. 3). Note that simulations of various geometrical cloud 
configurations can be constructed so that Ng can be empirically derived as a function of N, 
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b. Fluxes 


A consequence of the finite nature of each cloud in the array is to increase the effective cloud 
fraction because N, is always greater than N. This implies that energy loss to space from a broken 
cloud array is less than that computed using tiie cloud fraction without considering tlie sides of tiie 
clouds. 

The effective black body emission temperature at 10 iim from a region containing the model 
cloud anay at 2SS*K overlying a black ground at 290°K is shown in Fig, 5. The curve marked a ■ 0 
is the planiform case. It is evident that errors in flux of 10'^' equivalent temperature can be made 
for clouds of unit aspect ratio. This significantly affects the radiation balance because radiative models 
incorrectly weight clear and overcast fluxes as in (11). If a significant portion of the global cloud cover 

is broken, then (*^out)cuboidal ^ ^planiform’ ^ especially important for simulations of the 
Earth Radiation Budget Experiment. 

c. Net Flux Divergence 

Of possible impact in dynamic models is the fact that the cooiing rate across the cloud layer is 
different from the planiform case because of the radiation intercepted by the sides of tlie cloud 
(EUingson and Kolczynski, 1980). If a regular array of black isothermal clouds at Tq is embedded 
in a non*participating atmosphere over a ground at T^ , then the difference in net flux between the 
top and bottom of the cloud layer is 

AF„«-ir8(2Bo-Bi) (13) 

where 5 is the effective emittance of the layer and B^ and Bj are the Planck functions at T^ and Tj 
respectively. It follows from (13) that the cooling rate across the layer, calculated for the cuboidal 
cloud model is related to that obtained using the planiform assumption by 

^^^net^cuboidal _ ^ 

^^^net^pUniform ^ 
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Tliis quantity is plotted in Fig, 6 for various aspect ratios as a function of th^ cloud fraction, N. It 
is evident that for cumulus cloud fields with N » 0.3, which is a typical value, the cooling rate in the 
atmospheric window could be two to three times that derived from the plaiuform cloud assumption. 

d. Radiances 

If the array is composed of black isothermal cubic clouds, the radiance in any direction depends 
on the field of view obscured by thfe clouds. When viewed from above, only the cloud frucilon N ob- 
scures the ground, whereas from any lower observation angle, the sides of the cloud also obscure the 
ground. This effect is sliown in Fig. 7 for two viewing directions*, ^ ■ 0* , in which the observer views 
parallel to the array seeing only the top and one side face of each cloud, and ^ ■ 4S”, which is a diag- 
onal view in which the top and two adjacent side faces are seen. Individual curves are marked with 
the cloud fraction; the aspect ratio is ur;illy. The cuboids are at 2S5°K while the ground is at 290”K. 

The radiance profile is governed exclusively by geometry; the sharp kink in the individual curves 
is at the angle beyond which the underlying ground is completely obscured in the viewing direction. 
At higher zenith angles, for ^ » 0”, the radiance is 

I-\/NBo+( 1 -VbTjBi (15) 

Only two cases are shown for 0 ■ 45” because the geometrical arrangement is very complex, but that 
is sufficient to illustrate the effect of viewing direction. The radiance drops off much more rapidly 
for this case. Figure 7 suggests that the brightness temperature measured by a satellite instrument 
cannot be related to the cloud top temperature for broken cloud fields without considering the 
geometry of the array and the viewing direction if the field of view includes a number of cloud 
elements. 

3. Monochromatic Radiative Transfer through an Array of Non-Black Gouds 
a. Isothermal Gouds 

Let us now consider a regular array as in Fig. 1 but we introduce the requirement that the cloud 
is a cuboid of water droplets with a size distribution characterized as C-I by Deirmendjian (1969). 
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The radiative properties at 10 /urn of such cloud droplets are computed from Mie theory. Unlike the 
black blocks of Section 1, these clouds are partially transparent to 10 /um radiation. It is also assumed 
that the .'louds in the array are identical and isothermal while the ground beiow is non-retlecting. It 
will be shown that the results of the previous section can be used to obtain fluxes and radiances 
emanating from such an array, 

The thermal radiation emerging from cuboidal clouds was determined by means of tire VATS 
technique of Harshvardhan et al. (1981). It is sufficient to note here that fluxes leaving the faces of 
the cuboid may be expressed in terms of the mean radiance, Ig * 1 /4ir Id£2, where SI is the solid 
angle. Ig is obtained from the solution of 

v2|g-\2(i^.Bj (16) 

where 

\-kI3(l -w)(l 

with k, u and g being the extinction coefficient, single scattering albedo and asymmetry parameter, 

•• 

respectively. 

b. Boundary Conditions 

The boundary conditions for the problem are defined by the fluxes incident on each face. In 
general, the flux in any direction is a function of position, say F(x, y, z) and in particular, the flux 
exiting each side face is a function of position on the face. For a cloud of equal sides, s, and height 
2 *, with the origin at the center of the base, we have (Pxf,x+K»t/2 “ nomenclature 

used here is that of Davies (1978) and Harshvardhan et al. (1981). The first subscript identifies the 
face while the second subscript identifies the direction. 

The flux boundary conditions at the top and bottom faces are as for the single cloud, i.e. 

^^z+,z-^z*z* * ® 
and 

(Pz-,z+)z-o * ^^1 
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with B| the ground emission. The average tlux exiting tlie .side face at x s/2 is 


(^x+,x+)x-s/2*'~* ^ ^>9) 

If the regular array is made up of identical cuboids, then by symmetry, * (^t,y±>y-±s/2. 

and this may now be used to specify the side boundary conditions. 


TIte side faces receive radiant flux from neighboring clouds as well as that portion of tlie ground 
that is not obstructed by those clouds. If we assume that the flux exiting the side faces is isjtroplc 
and uniform over the face, it is possible to compute the incident flux on each side face of the cloud. 
Let A| be the area of each side face and A 2 the total area of all the faces visible from this face for a 
particular cioud fraction. Figure 2 shows that A 2 varies with the configuration because of partial ob> 
scuration by nearby clouds. The energy per unit area falling on Aj that is diffusely emitted by A 2 
^ ^^x+,x+^x*'$/2 ^21 “"8^® factor between Aj and Aj defined as the 

fraction of energy diffusely emitted by A 2 falling on Aj . By reciprocity (Sparrow and Cess, 1978) 


®21 ^2 


'12 


( 20 ) 


and 3 12 is identical with 3 of Eq. (4) et seq. Therefore, to the extent that the cloud faces may be 
approximated by isotropic emitters of uniform flux, the side boundary condition for the array prob* 
iem is written as 


CF, 


ttB, 


x+,x-^x*»/2 “ ^(i^x-,x-^x »-»/2 


( 21 ) 


In (21), (Fx+^x-^x*s/ 2 “ ^® inconiing flux at the x * s/2 face whereas (Fx_^x-^x»-s/2 '^®“" 

exiting flux out of the opposite face at x » -s/2. 


It is necessary to examine the two assumptions made in deriving (21). The assumption of iso- 
tropy can be justified to a certain extent if the cloud is sufficiently thick, kz* > 10, so that the up- 
ward and downward components of the mean side flux are nearly equal. This is because the radiant 
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flux is emission dominated. Moreover, the problem would become far more intractable if the angu- 
lar distribution were to be accurately considered. 

The assumption of uniform flux over the side face is not correct for single clouds (see Harshvard- 
han et al., 1981). As tlie cloud fraction increases and the spacing decreases, an error will be introduced 
by assuming that the mean flux is incident on the neighboring cloud. In the limit that the cloud sides 
touch each other, the input to one side face is the output of the opposite face, point for point, i.e. 

^^x+,x-\»s/2 “ ^^x-,x-^x"-s/2 (22) 

Note that if we set 3 « 1 in (21 ), the R.H.S. of the equation contains the mean flux and not 
the flux at each point. 

Figure 8 illustrates the exchange geometry. We merge (21) and (22) such that it a applicable 
over the entire range of cloud fractions. The weighted boundary condition is 

_ irB, 

^^x+,x-^x»*/2 * ^^x-,x-^x»-«/2 ’^^^x-,x-^x*-s/2^ 3 (1 " 2 ’ 


^ is chosen to weight the incident flux; it is the ratio of two angle factors, 

^Aj-dAi j2 




^dAj-dAj 


; 4rrG(a,N) 


where the funct'on 


G(a,N)^ 


J (a/2R) 


Wr+(a/2R)“ 


tan 


-1 


(1/2R) 


LVl +(a/2R)^J 


(1/2R) . 

+ - — tan * 


Vl +(l/2Rf 


(a/2R) 


yi +(1/2R)^J 


(24) 


(25) 


with R * d/s * 1/VN “ I and a ■ z*/s, the aspect ratio (Tltis expression is given in the Appendix of 
Sparrow and Cess, 1978). An infinitesimal area dAj receives flux F from an element dAj exactly op- 
posite it and flux F from the whole opposite face Aj. Figure 9 shows that 1 // tends to unity as N be- 
comes small and tends to zero as the faces approach each other. The precise weight in boundary con 
dition (23) and choice of \jj is somewhat arbitrary but is physically reasonable. The contribution to 
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the flux arriving at a point on the sides of the clouds from any point on the opposite faci,<^ ^'iepends on 
the angle factor between the points. In an approximate sense, the mean flux may thus be identified 
with the angle factor between the opposite face as a whole and the receiving point. Tlie angle factor 
^A-dA- on *he location of dAj but the center Of the face has been used to position dAj. 

For the special case w.hich exliibits uniform flux emission, F * P, and (23) reduces to (2 1 ). 

c. Fluxes 

The total flux emanating from the array and underlying ground can now be obtained using the 
single cloud solution obtained with the boundary condition (17), (18) and (23). If (^zi-.z-t-^z-z* 
the mean flux out of the top face of each cloud, then the power emitted is A^fF^^ and if 

^*^x+,z+^x»s/2 ^ upward component of the mean flux out of each side face, then the power 
emitted to space by ail four sides is 4Aj(Fjj^ Therefore, the upward flux averaged 

over this array may be written as 

Fou, - N<P„„, H-4afl- 3) 4 (1 - N.) ,rB, (26) 

If the clouds are considered black, (^x+,z+^x«s/2 “ ^/^^^x+,x+V»s/2 “ (26) reduces to (10). 

I 

Computations have been made for a » 1 (cubic cloud) and the radiative model mentioned earlier 
for which the single scattering albedo, » 0.638 and asymmetry parameter, g * 0.865 at 10 jam. Fig- 
ure 10 shows the average upward flux from a regular array of cubes of optical thicknesses on each 
side varying from 0.5 to 20, Results are presented as an effective black body temperature. Tliese 
results may be compared with the curve marked a * I in Fig. 5 to note the departure from the black 
case. It is evident that an optical depth of 10 or greater may be approximated quite well by black 
cuboids. 

d. Radiances 

Figure 1 1 shows the radiance at ■ 0® for various cloud fractions for a regular array of cubic 
clouds of optical dimensions (ks, ks, kz*) » (10, 10, 10). Tire ground and cloud temperatures are 
the same as before and the plot may be compared with Fig. 7 which is the black case. It should be 
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stressed here that these solutions are not exact but VATS approximc.tions. However, Harshvardhan 
et al. (1981) have shown that the approximation is good, especially for large optical depths and these 
non*black cuboidal cloud results appear to be reliable. Tlie black assumption is very good for zenith 
angles less than 60° as can be seen from Fig. 1 2 wiuch shows the difference between the brightness 
temperature fc ' black and non-black clouds with a * 1 . The maximum error in radiance is less than 
2°K when 6 < 60° for all arrays of optically tliick clouds. Tliis situation is frequently encountered 
in remote sensing problems. 

The simulated radiance from partly cloudy fields is of the form 

I-IpN + B, (1 -N) (27) 

where Ip is the planiform Elution, Bjthe ground contribution or clear sky solution and N represents 

the cloud fraction. If N is taken to be the cloud fraction when viewed normally, (27) will be erron- 
eous for non-planiform clouds. To illustrate this, we liave computed I according to (27) and com- 
pared it with the present solution for an array of cubic clouds of optical dimensions (10, 10, 10). 

The value of Ip used in (27) is the corresponding planiform solution for optical depth 10. The dif- 
ference is plotted in Fig. 13. The cloud is isothermal in both cases at Tg * 2SS°K while the ground 
temperature is Tj ■ 290°K. It is evident that (27) is not valid for any zenith angle and cloud fraction 
with errors of 4°K in brightness temperature for a typical field of N = 0.3 viewed at 0 * 30°. For 
fair weather cumulus, with a ~ 1 , the geometry is more significant than the detals of the radiative 
transfer microphysics and good results can be obtained with the black block assumption. 

e. Non-Isothetmal Clouds 

We have only considered isothermal clouds so far; however, it can be demonstrated that iso- 
thermality is not a serious drawback for fair weather cumulus and other cloud formations that are 
not extensive in the vertical direction. Consider a cuboidal cloud array with a » 1 and optical dimen- 
sions (10, 10, 10). A ver^iical optical depth of 10 at 10 Atm corresponds to a maximum vertical dimen- 
sion of about 1 km (Deirmendjian, 1 969), so the temperature difference between cloud top and base 
would be about 5°K. We therefore compute the radiance at <> = 0° from a non-black array of cuboids 
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(c5* 0.638, g» 0.865) with cloud top temperature, » 255“K, cloud base temperature, ■ 260®K, 
overlying the ground at Tj ■ 290° K and compare this with the isothermal black cuboid radiance with 
To - 290°K. 

The difference curves are shown in Fig. 14, which is very similar to Fig. 1 2, the isothermal com- 
parison. Again, the maximum error is only a little more than 2°K for B < 60°. One reason that the 
lapse rate is not very important for sufflciently opaque clouds is that for small values of d, the bulk 
of the radiation is emitted from the top and for larger zenith angles, only the upper portion of the 
sides are visible, the rest being obscured by neighboring clouds. !n so far as clouds are optically thick, 
replacing a broken cloud array by black emitters at the same temperature as the cloud top is superior 
to the use of (27) in computing exiting radiances. 

4. Summary 

The results presented here show that it is possible to model the infrared radiation emanating 
from a broken cloud field if only the geometry of the array and the cloud top temperature are de- 
fined. Although a very simple regular array of cuooids was considered, this study compares a 
black cloud and a scattering-absorbing water cloud with appropriate radiative properties. The non- 
black problem was solved by extending the single cloud solution as obtained by a two-stream 
approximation to an array solution with mutual interaction. One conclusion of this study is that 
clouds may be approximated by black opaque finite emitten in the infrared if the optical dimensions 
exceed 10 on each side. This will be particularly useful in modeling irregular arrays of non-uniform 
sized clouds. It was also shown that a simple area-weighting of overcast and clear sky radiances for 
broken cloud field radiative simulations is highly erroneous when the aspect ratio of the cloud ele- 
ments is about 1 or larger. 
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Figure 2. Top view of array showing cloud faces visible 
from hatched face on element 00. 
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Figure 3. Effective cloud fraction for black bar clouds. Comparison 
of theory and measurements is for a > 1 . 
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CLOUD FRACTION, N 

Figure 4. Effective cloud fraction for black cuboidal clouds, Ng, as a 
function of cloud fraction, N, for aspect ratios marked on curves. Full 
lines are given by Eq. 1 2 in text. 





CLOUD FRACTION, N 

Figure 5. Upward flux from an array of black cuboids expressed as effective 
black body temperature for Tq » 255®K and Tj » 290®K. Aspect ratios are 
marked on individual curves. 





Figure 6. Ratio of cuboidal ai 
as a function of cloud fractior 
ual curves. 
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Figure 8. Geometry for energy exchange between opposite 
faces showing interchange between a point and the entire 
opposite face, and the opposite point. 



CLOUD FRACTION, N 

Figure 9. The function defined by Eq. 24, \p vs N, the cloud fraction. 
Aspect ratios are marked on individual curves. 
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Figure 1 0. Same as Fig. 5 but for non-black clouds of unit aspect ratio. 
Radiative properties are w » 0.638, g = 0.865. Optical thickness on each 
side of the cuboids is marked on individual curves. 
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Figure 1 1 . Same as Fig. 7 but for non-black clouds of unit aspect ratio 
at 0 * 0°. Optical dimensions are (10, 10, 10). Radiative properties are 
w - 0.638, g * 0.865, Tq • 255”K, T, » 290“K. 
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Figure 1 2. Difference in brightness temperature at 0 * 0° between an 
array of isothermal black cuboids of unit aspect ratio and non-black 
cuboids of optical dimensions (10, 10, 10). Radiative properties of 
non-black cloud are <2 = 0.638, g = 0.865. Cloud fractions are marked 
on individual curves. Tq * 255"K, Tj = 290“K. 





Figure 13. Difference in brightness temperature between a planiform 
cloud of optical thickness, kz*, equal to 1 0 and an array of isothermal 
non-black cuboids of optical dimensions (10, 10, 10) at ^ = 0". Radi- 
ative properties for both clouds are w = 0.638, g = 0.865. Cloud frac- 
tions are marked on individual curves. ® 255®K, Tj * 290® K. 
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Figure 14. Same as Fig. 12 but with cloud top, Tq. = 255'’K; 
cloud base, = 260“K; ground, Tj = 290^K. 
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